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Three dimensional topological insulators embody a newly discovered state of 
matter characterized by conducting spin-momentum locked surface states that 
span the bulk band gap as demonstrated via spin-resolved ARPES measure- 
ments [l-8|. This highly unusual surface environment provides a rich ground 
for the discovery of novel physical phenomena j5l-l26l| . Here we present the first 
controlled study of the topological insulator surfaces under strong Coulomb, 
magnetic and disorder perturbations. We have used interaction of iron, with a 
large Coulomb state and significant magnetic moment as a probe to systemati- 
cally test the robustness of the topological surface states of the model topologi- 
cal insulator Bi2Se3. We observe that strong perturbation leads to the creation 
of odd multiples of Dirac fermions and that magnetic interactions break time 
reversal symmetry in the presence of band hybridization. We also present a the- 
oretical model to account for the altered surface of Bi2Se3. Taken collectively, 
these results are a critical guide in manipulating topological surfaces for probing 
fundamental physics or developing device applications. 

Bismuth selenide has been discovered by angle resolved photoemission spectroscopy 
(ARPES) to be a topological insulator with a bulk band gap of about ~300 meV 7H9|. 
Spin resolved photoemission studies reveal that surface electrons in Bi2Se3 form a Dirac 
cone spanning the bulk insulating gap, composed of spin-momentum locked helical states 
revealing its topological order (Fig. la). The Fermi level of grown crystals is usually found 
to be located in the bulk conduction band due to selenium vacancy defects, however it was 
subsequently shown that with Ca doping or NO2 surface deposition, the Fermi level can be 
placed at the Dirac point reaching the topologica 



interactions can be controlled via Fe or Mn 



transport regime [7|, [lOj] and magnetic 
Placing the Fermi level at the Dirac 



point in the presence of magnetic impurities in the bulk can lead to a small gap at the Fermi 
level (Fig. lb), however the full character of this gap cannot be decisively resolved with 
the current resolution of ARPES due to lineshape broadening effects. So far, no systematic 
surface deposition of magnetic impurities to elicit a large systematic magnetic response and 
bring about controlled changes in the surface band structure has been explored. Probing 
the effect of magnetic perturbation on the surface is more relevant for potential applications 
than the previous studies of bulk dopants, because topological insulators need to be in con- 



tact with large moment ferromagnets and superconductors for device applications 13l-ll9|. 
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The focus of this Letter is the exploration of topological insulator surface electron dynamics 
in the presence of magnetic, charge and disorder perturbations from deposited iron on the 
surface. 

Bismuth selenide cleaves on the (111) selenium surface plane, providing a homogeneous 
environment for deposited atoms. Iron deposited on a Se^" surface is expected to form a 
mild chemical bond, occupying a large ionization state between 2"*" and with roughly 4 



fiB magnetic moment 



27j . As seen in Fig. Id the surface electronic structure after heavy 



iron deposition is greatly altered, showing that a significant change in the surface electronic 
environment has been achieved. Five (odd number of) surface bands intersect the Fermi level 
rather than just one, and extend below the Fermi level in the form of multiple Dirac cones. 



Unlike the bulk electronic states of Bi2Se3 



24j |. these features have no z-axis momentum 



dependence, confirming that they are largely two dimensional in nature (Fig. le). 

Figure 2 shows the surface evolution as a function of deposition time in crystals with 
chemical compositions of (Sample #1) Bii.9975Cao.oo25Se3 and (Sample #2) Bi2Se3 varied to 
tune the carrier density. The rate of iron deposition was similar to 4% of a monolayer per 
minute, as detailed in the SI, and deposition data for an additional sample with intermediate 
carrier concentration (Sample #3) are also presented in the SI. The bulk electronic state of 
Sample 7^1 prior to deposition is slightly p-type with significant resistivity {p = 23 mficm), 
and Sample #2 is electron doped from Se vacancy defects in as-grown Bi2Se3. In each 
case, it is observed that the presence of positively charged Fe surface ions progressively 
lowers the energy of the surface state and causes the appearance of new surface states with 
energy/momentum contours similar to the bottom of the bulk conduction band. 

The changes brought on by iron deposition can be seen most strikingly in measurements 
on the bulk p-type Sample #1. After approximately six minutes of deposition, new surface 
states that do not exist in unperturbed samples become visible at the Fermi level, with new 
Dirac points labeled Dl and D2 above the Dirac cone that is seen in as-grown Bi2Se3 (DO). It 
is also after approximately six minutes of deposition that a gap begins to be apparent at the 
DO Dirac node, evidenced by a parabolic shape near the Dirac node, separating the upper 
and lower Dirac cones of the original surface state. This gap can be seen clearly at different 
incident photon energies in Fig. 2d, confirming that it is a feature of the surface and not the 
bulk electronic band structure. Electron velocities (band slope) near the Dl and D2 Dirac 
points increase monotonically as iron is added, showing that iron is increasing the "Rashba" 
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interaction term {{kx z) -a, with a representing the Pauh matrices) identified in theoretical 
models jsj. The number of surface bands intersecting the Fermi level between the F- and 
M-points progresses from one to three to five, with one band contributed by the original 
(DO) Dirac cone and two more bands contributed by each of the new (D1,D2) Dirac points. 
This is consistent with the Mod(2) character of surface electrons on a crystal with bulk 
topological insulator order, that the topological surface likes to maintain an odd number 
of Dirac states. After 12 minutes of deposition, the binding energy of the DO Dirac point 
was found to have sunk approximately 0.6eV in energy, and the electron binding energies 
ceased to change under additional deposition. When the chemical potential is positioned 
above the bulk conduction band minimum, as in Sample 7^2 (Fig. 2c), the dispersion of new 
surface states across the full bulk band gap is visible within the photoemission image. A 
new, strongly split surface band is observed in Sample #2 after five minutes of Fe-deposition 
with a (Dl) Dirac-point at the F-point, but no further (e.g. D2) bands appeared after longer 
deposition. 

Theoretical simulation of non-magnetic surface Coulomb perturbation on Bi2Se3 is shown 
in Fig. 3a, and qualitatively reproduces the progressive appearance of new Dirac points with 
increasing iron deposition. Through comparison with our numerical result, we can see that 
the experimentally observed surface states begin to pair off at momentum separation greater 
than ~0.lA~^ from the Brillouin zone center, with the upper DO Dirac cone approaching 
degeneracy with the lower Dl band, and the upper Dl band connecting to the lower D2 
band. The partner-swapping connectivity observed in the simulation and data is a simple 
way by which new states can be added to the surface band structure without disrupting 
the surface conditions required by the bulk topological insulator order of Bi2Se3 jsl. The 
spin-splitting of topological surface bands is often discussed as a special case of the Rashba 



effect (e.g. Ref. [3|, l22|), in which surface electronic states become spin-split by an energy 



proportional to their momentum k. Our data and simulations show that this description 
is only accurate for Bi2Se3 in a small part of the Brillouin zone surrounding the Brillouin 
zone center, because at momenta further from the F-point the electronic states pair off and 
are nearly spin degenerate (see Fig. 3b). This can be understood because the origin of 
;he topological insulator state in Bi2Se3 is a symmetry inversion that occurs at the F-point 
sl, H], and the electronic states close to the Brillouin zone boundary are similar to those of 
topologically trivial materials. The role of magnetic domains in reshaping low energy band 
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structure is expected to be more subtle than Coulomb perturbations, and may be limited 
to changes near the Dirac points and in self energy lineshape effects (see Fig. 3c-d). 

Recent theoretical studies suggest that the physical environment of magnetic impurities 
on a topological surface is very different from the surface environment provided by a nor- 
mal semiconductor such as silicon When a non- magnetic crystal is doped with 
magnetic impurities, long-range magnetic ordering can come about as a result of itinerant 
electrons exchange-mediating the magnetic interaction. In a normal three dimensional ma- 
terial, electrons that interact with surface-deposited magnetic impurities are free to scatter 
away from the surface over a 2tt solid angle, and the two dimensional magnetic interactions 
therefore typically decay over several Angstroms (e.g. Ref. 29(|). In topological insulators, 
surface state electrons are naturally confined to the surface in two dimensional Dirac cones, 
and it is suggested that interactions between deposited impurities can be mediated over 



many nanometers (see illustration in Fig. 4c) 



13 



28| . For broader theoretical investigations 



and experiments with nanoscale engineering capabilities, it is interesting to note that exotic 
magnetic phases (e.g. helical) may occur within certain specific configurations of magnetic 
impurities such as isolated ID chains, or from instabilities that arise as the chemical poten- 



tial is moved far from the surface Dirac point 



22 



26| . In the present case of homogeneous 2D 



deposition with chemical potential near the DO Dirac point, the momentum-locked spin po- 
larization of topological surface electrons is thought to support ferromagnetic order with an 
out-of-plane bias [B along the ±i direction) at high impurity densities 13|, llJ], unlike direct 
dipole-dipole interactions which are much weaker and favor in-plane magnetic orientation. 

Our numerical simulations have shown that the new Dl and D2 Dirac point electrons are 
localized deeper inside of the material than DO, and are expected to interact more weakly 
with electronic orbitals of the surface-deposited iron. Therefore, we focus on the DO electrons 
to understand the effect of iron magnetism on the topological surface. A gap appears at 
the DO Dirac point after heavy iron deposition, and the lower DO Dirac cone acquires a 
buckled shape with a local energy minimum at zero momentum. The mass induced in the 
upper DO Dirac cone after full deposition is approximately 0.1 Me (electron masses), and the 
surface band gap is similar to 100 meV. It is difficult to identify the exact time at which the 
jand gap appears, likely because there is no long-range order of iron spins on the surface 



30 



3l|, and differently ordered domains will yield differently gapped contributions to the 



photoemission signal. 
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Global ferromagnetism breaks time reversal symmetry, making it possible to induce a gap 
at the Dirac point, which is otherwise disallowed by the crystal symmetry. Figure 4a shows 
the result of modifying a Generalized Gradient Approximation (GGA) numerical prediction 
of the DO surface state by adding perturbative coupling to a surface layer of magnetic 
impurities with ferromagnetic out-of-plane order (equivalent to the Zeeman effect), yielding 
a dispersion that closely matches the experimental data. Modifying the surface further 
with negative-ion deposition could place the Fermi level inside of the surface band gap (Fig. 
4a(right)). Details of this calculation and factors that can modify the surface state dispersion 
are discussed in the SI. If the DO band gap were due to an isotropic out-of-plane magnetic 
field, a similar gap would be expected at the Dl and D2 Dirac points, which is clearly not 
observed. Because magnetic field perturbations from the sub-monolayer iron coating are 
quite weak, the magnetic phase term in perturbative Hamiltonians can be discounted, and 
the dominant magnetic effects are expected to be described by a local Heisenberg exchange 
term representing direct wavefunction overlap (hopping) between surface electrons and spin 
polarized Fe 3d orbitals 13|. Evidence of Landau levels is neither seen in our data nor in 
this context expected from theory. The general effect of Zeeman-type magnetic symmetry 
breaking on topological Dirac surface electrons is shown in Fig. 3c, illustrating that out- 
of-plane magnetic order will cause a gap, and in-plane magnetic domains will cause bands 
to shift in momentum space jl^. Based on the appearance of a gap at DO in our data, we 
conclude that much of the iron-deposited crystal surface is occupied with domains that have 
net out-of-plane magnetic moment. 

Contrary to the usual trend of deposition experiments, in which photoemission images 
become increasingly blurry as molecules are haphazardly added to the surface, we identify a 
regime in which the image becomes qualitatively sharper with the increasing coverage of iron 
(seen in in Fig. 4d and panels 1-3 of Fig. 2a), corresponding approximately to the increasing 
clarity of the DO gap. Momentum-axis width of bands {6k) was self-consistently measured 
through Lorentzian fitting at binding energies between the DO Dirac point and the onset of 
the Dl band structure. Momentum width is inversely related to the electronic mean free 
path, and the reduction in 6k as the gap appears is likely indicative of a magnetic disorder-to- 
order transition, leading to a reduction in scattering. Reducing the prevalence of magnetic 
domains with in-plane moment could also cause the surface state bands to appear narrower 
even if electronic mean free paths are unchanged, as can be observed through comparing the 
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domain- averaged band profiles in Fig. 3d. The signature of large ferromagnetic domains with 
in-plane moment is a twin-peak feature that is not seen in the data. The close correspondence 
between reduced momentum-distribution-curve (MDC) or 5k width and gapping of the DO 
Dirac point therefore strongly indicates the occurrence of an out-of-plane magnetic ordering 
phase transition after approximately six minutes of iron deposition. 

This apparent ordering transition driven by magnetic interactions that are mediated by 
the topological surface state is an indication of how strong topological insulator order changes 
the physical environment of the material surface. The observations reported here, includ- 
ing out-of-plane surface magnetism and the appearance of new Dirac surface states, open 
a window into how topological surface states are formed and interact with various pertur- 
bations. As such, they are significant for theoretical understanding of the formation of the 
topological insulator state in particular materials, and have direct implications for proposed 



devices utilizing magnetic or Cou 



as magneto-electric junctions 14 
topological transistors. 



omb-charged interfaces with topological insulators, such 



17j and capacitors 



18] to integrate into multifunctional 



Methods summary: Angle resolved photoemission spectroscopy (ARPES) measure- 
ments were performed at the Advanced Light Source beamlines 10 and 12 using 35.5-48 
eV photons with better than 15 meV energy resolution and overall angular resolution bet- 
ter than 1% of the Brillouin zone (BZ). Samples were cleaved and measured at 15°K, in a 
vacuum maintained below 8x10^^^ Torr. Momentum along the z axis is determined using 
an inner potential of 9.5 eV, consistent with previous photoemission investigations of un- 
doped Bi2Se3 jsl, Q . Fe atoms were deposited using an e-beam heated evaporator at a rate 
of approximately O.lA/minute. A quartz micro-balance supplied by Leybold-Inficon with 
sub- Angstrom sensitivity was used to calibrate the iron deposition flow rate. Adsorption of 
NO2 molecules on Caa;Bi2_a;Se3 was achieved by controlled in situ exposures under static 
flow mode, with care to minimize photon exposure of the adsorbed surface. Large single 
crystals of Ca-x^h-x^es were grown using methods described in Ref. [10|. Surface and bulk 
state band calculations were performed for comparison with the experimental data, using 
the LAPW method implemented in the WIEN2K package 32|. Details of the calculation 
are identical to those described in Ref. [8j. 
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FIG. 1: Iron deposition strongly modifies the topological surface: a, Uniformly electron 
doped Bi2Se3 has a single surface state Dirac cone, b, When the surface chemical potential of 
as-grown Bi2Se3 is lowered to the Dirac point by NO2 deposition, we observe a slight gap in the 
leading edge of ARPES intensity, c, The hexagonal surface Brillouin zone of Bi2Se3 is drawn above 
a diagram of the three dimensional bulk Brillouin zone, d, A second derivative image of new surface 
states in Bi2Se3 (Sample #1) after surface iron deposition is labeled with numerically predicted 
spin texture from Fig. 3b. e, Low energy features from (d) have no z-axis momentum dispersion 
as seen from the data taken with varying incident photon energy (37 eV to 29 eV) confirming the 
two dimensional character of the state. 
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FIG. 2: Iron doping creates five Dirac cones: a, Progressive surface doping of optimally 
insulating Ca2;Bi2_xSe3 (Sample #1) causes the successive appearance of two new Dirac points 
("Dl" and "D2"). b, Second derivative images of the Fermi surface and low energy surface states 
are shown for Fe-deposited Sample i^l. c, Iron deposition on as-grown n-type Bi2Se3 (Sample #2) 
results in a similar final spectrum to panel (a), but with only a single new Dirac point, d, The 
iron- induced DO gap is observed at different binding energies for Sample ^^1 and (inset) in a second 
derivative image. Existence of the surface band-gap is not sample dependent. 
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FIG. 3: Surface state simulations: a, A first principles GGA model of Bi2Se3 surface band 
structure (at left) exhibits new Dirac points when orbital energies are lowered by 0.68eV on (center) 
the first and (right) the first two quintuple layers of the lattice. Surface state energy changes 
nonlinearly with this perturbation. Inset schematics show the Fermi surface evolution (not to 
scale), b, Bands from panel (a, right) are labeled with t for right handed and J, for left handed 
chiral spin texture, c, The effect of Zeeman coupling on a surface state Dirac cone defined by a 
Rashba Hamiltonian is considered for different magnetic field orientations, based on the surface 



state Hamiltonian in Ref. [141 ]. Red arrows illustrate the path of a constant energy ARPES 
measurement, corresponding to profiles in panel (d). d. The domain-averaged momentum axis 
distribution of a single surface state band is shown for a surface populated with (black) out-of-plane 
ferromagnetic domains, (red) in-plane domains and (blue) domains with arbitrary 3D orientation, 
assuming a band velocity of 3.4eV-A and magnetic interactions strong enough to open a O.leV gap 



when oriented along the z-axis. 
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FIG. 4: Surface evolution under disorder and magnetism: a, (left) A z-axis magnetic 
perturbation causes (green) GGA predicted surface states to conform to tlie iron doped dispersion, 
(right) A surface insulator phase can be obtained by adding ~4 Langmuirs [7] of non-magnetic 
NO2 on the surface subsequent to Fe deposition, b, A cartoon illustrates Fe deposition on the Tl 
surface, c, (black arrows) Electron trajectories in TIs adhere to the surface, helping to mediate 
2D magnetic interactions. Red arrows represent Fe spin, d, The development of a gap at the 
lowest energy ("DO") Dirac point is compared with (bottom) the half maximum peak width (6k) 
of connected states in the DO upper Dirac cone. Error bars represent estimated uncertainty based 
on the results of different fitting techniques. 
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